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Abstract. We reinterpret and generalize conjectures of Lam and Williams 
as statements about the stationary distribution of a multispecies exclu- 
sion process on the ring. The central objects in our study are the mul- 
tiline queues of Ferrari and Martin. We make some progress on some of 
the conjectures in different directions. First, we prove their conjectures 
in two special cases by generalizing the rates of the Ferrari-Martin tran- 
sitions. Secondly, we define a new process on multiline queues, which 
have a certain minimality property. This gives another proof for one of 
the special cases; namely arbitrary jump rates for three species. 



1. Introduction 

We study a totally asymmetric simple exclusion process (TASEP) on the 
ring Z/NZ. There are particles of n different classes present and every 
position of the ring is occupied by exactly one particle. It is a continuous 
time Markov chain in which a particle can jump over (i.e. trade place with) 
a particle from a larger class to its left. Lam [9] conjectured properties of 
this Markov chain which he needed for his work on infinite reduced words of 
affine Weyl groups. This TASEP has been studied independently by several 
authors and in particular it has been given a beautiful solution in terms of 
so called multiline queues by Ferrari and Martin [8J, which we describe in 
Section [2j These intricate objects give a solution to, and predate, one of 
Lam's conjectures about the partition function, but more work is needed to 
resolve others, in particular the stationary weight for the identity. 

In [TU] Lam and Williams generalized the model, but in a different lan- 
guage, by including different jump rates for different classes of particles, 
so particles of class i jump with rate X{. They obtained surprisingly nice 
stationary weights, which they conjectured to be polynomials in the jump 
rates with positive integer coefficients, see Conjecture |2.5| The purpose of 
the present paper is to report some advances on these conjectures. 

In particular, see Section [3j we solve this inhomogeneous TASEP with at 
most 3 different kinds of particles, Theorem |3.2[ verifying this conjecture. In 
fact we do this in two different ways, first using the transitions for multiline 
queues by Ferrari and Martin and then by defining a new different Markov 
chain, see Section [5| which also projects down and gives a solution of the 
TASEP with three different classes of particles. We include this second proof 
because we have higher hopes that it may generalize to n > 4. We also prove 
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in Section [4] that the Markov chain defined by Ferrari and Martin is enough 
to understand the power of x\, setting xi = • • • = a%_i, for an arbitrary 
system with exactly one first class particle. The interested reader can verify 
the statements in the paper by downloading the Maple package InhomTasep 
from either the homepage of one of authors or the arXiv source. 



1.1. Background on TASEP. The general exclusion process can be de- 
fined on an arbitrary graph. One starts with a configuration of particles on 
the vertices of the graph, where every vertex can be occupied by at most 
one particle. The process involved hopping of the particles according to a 
Markovian rule, and the exclusion condition refers to the fact that at most 
one particle can be at any site. 

The exclusion process was probably first studied in the biological litera- 
ture. It was proposed as a prototype to describe the dynamics of ribosomes 
along RNA |13j . They were studied systematically by probabilists in the 70s 
starting with the work of Spitzer [14], who coined the term. A lot of results 
are now known about the exclusion process on Z [12J. 

Two starting points for a more combinatorial understanding of exclu- 
sion processes on a finite state space was for the TASEP |5J by Duchi and 
Schaeffer and for the PASEP j2] by Corteel and Williams. There are many 
subsequent papers relating these to well-known combinatorial constructs. 
Recently, this approach has led to the first combinatorial interpretation of 
the moments of the Askey- Wilson polynomials [3J. 

Motivated by questions in statistical physics, Derrida, Janowsky, Lebowitz 
and Speer [3] obtained the stationary distribution for the TASEP on Z/iVZ 
with two species of particles, first-class and second-class, in addition to va- 
cancies. The general solution for the stationary distribution for any number 
of different classes of particles came from Ferrari and Martin [8], who built 
on their own work on multiline queues [7] and used ideas from Ferrari, Fontes 
and Kohayakawa [6 J and Angel jT]. 

The one- dimensional totally asymmetric exclusion process arose again, 
using very different terminology, in the study of random walks in Weyl 
alcoves by Lam [9] . He managed to prove that some results about the infinite 
Markov chain could be obtained by studying the stationary distribution 
of a finite one on permutations. This finite Markov chain turns out to 
be equivalent to the multispecies TASEP on a ring discussed above with 
one particle of each class. Further work along this direction led to more 
conjectures by Lam and Williams [10J, in particular the positivity conjecture 
addressed in this paper. 
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2. Multispecies Exclusion Processes and Multiline Queues 

The aim of this section is to explain the exact rules for the multispecies 
TASEP and for the multiline queues and then explain how the latter projects 
down to give a solution of the former. 

Before we define these processes, we quickly recall the part of the theory 
of Markov chains relevant to us. The continuous-time Markov chain on a 
(finite) state space Q is defined by the so-called transition matrix or M, 
whose rows and columns are labeled by elements of Q. For G\,C% 6 f2, we 
take the convention that the (C2, C\) entry is the rate of the transition from 
C\ — > C2 if Ci / C%. The diagonal entries (C\,Ci) is the negative of the 
sum of the rates of transitions leaving C±, not counting loops. 



This ensures that column sums are zero and consequently, zero is an eigen- 
value with row (left-) eigenvector being the all-ones vector. Assuming that 
the Markov chain is irreducible, it has a unique stationary distribu- 
tion jTT] . This is given by the entries of the corresponding column (right-) 
eigenvector v with eigenvalue zero. 

The stationary distribution is thus determined by solving the equation 
M.v = 0. This can be written in the following way. Let v(C) be the 
stationary distribution of configuration C, which can be viewed as the entry 
of the vector v at position C. Then 



This is the all-important master equation for C. We will use the unique- 
ness (up to scaling) of the solution of these equations to prove all our results 
about stationary distributions. 

2.1. Multitype TASEP on a ring. We will define the model on the ring 
Z/2VZ. The state space of the n-species exclusion process O m is defined by 
an n-tuple m = (mi, . . . , m n ) of positive integers where mi + • • • + m n = N. 
The configurations C G fi m are denoted by A-tuples with mi number of l's, 
?7i2 number of 2's, and so on, ending with m n number of ns. Clearly, 



Ci denotes the class of the particle at the ith site in the lattice. Each species 
can be thought of as a positive integer and lower integers have higher priority 
than higher ones. This might seem counterintuitive, but this notation makes 
sense if we think of first-class particles (of class 1) being higher in priority 
that second-class particles (of class 2). 

Definition 2.1. The multispecies exclusion process on Q m is defined 
by local transitions involving sites i and i + 1. If the current state is C , with 




te(C 2 ^C7i), 

rate(Ci ->• C), d = C 2 . 



rate(C7' -> C)v(C) = r ate(C -> C')v(C). 
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Ci = a, Ci+i = b, then 

I b a with rate 1 if a > b, 
I a o otherwise. 

In words, each particle carries an exponential clock which rings with rate 
1, and the particle tries to jump to its left whenever the clock rings. If the 
particle to the left is of less importance the jump takes place. Otherwise, 
the configuration is unchanged. 

Therefore l's always move left and n's always moves right. Particles whose 
class is between 1 and n sometimes move left, and sometimes right. We 
consider this to be a homogeneous model because the rate of all transitions 
is 1. The model of Lam and Williams |10j on permutations, on the other 
hand, is inhomogeneous because transition rates depend on the particles 
being interchanged. The careful reader will notice that our Markov chain is 
different from theirs because we interchange neighboring positions whereas 
they interchange neighboring values. In permutation group theory, this is 
just the difference between multiplying with simple transpositions on the 
right or on the left, so there is an easy bijection between the permutations 
in two chains. In particular, the set of stationary weights are given by the 
same expressions. 

We now generalize their model to multipermutations, where the rate of 
the exponential clock attached to particle of species j is Xj. They are to be 
thought of as positive numbers, but we treat them as formal parameters. 

Definition 2.2. The inhomogenous multispecies exclusion process 

on £l m is defined by local transitions involving sites i and i + 1. If the 
current state is C , with Ci = a, Ci+i = b, then 

, I b a with rate x^ if a > b, 
a b — > < (2) 
I a b otherwise. 

We give the m = (1,1,1) example in detail now. The graph of the Markov 
chain and the stationary weights of each configuration are given in Figure [TJ 
The transition matrix in the lexicographically ordered basis is given by 

\ 



X2 
Xl 

-Xl - X 2 J 

Notice that permutations which are obtained by rotations of one another 
have the same stationary weights. This is true in general because the tran- 
sition rates defined in Q depend only on the particle classes, not on the 
positions. 

2.2. Ferrari-Martin multiline queues. We now define the multiline queues 
of Ferrari and Martin [8j. As before, let m be the type of the process, with 
n species and N sites. Then the state space of the process, ^ s defined 

on a cylinder of circumference ./V and height n — 1, each site of which is 
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(xi+x 2 ) 231 




312 (X1+X2) 



(*i) 



Figure 1. Transitions of the Markov chain for n = 3. The 
stationary weights are given in parenthesis and in red. 



either occupied by a * or a o. We think of * as being occupied and o as 
being vacant. It will be convenient to set M% = Y2 l j=i m j- At row j, for 
1 < j < n, the configuration will contain Mj *'s and N — Mj o's. The rows 
are numbered from the top. At each row, all possible configurations with 
the prescribed number of *'s and o's are allowed. Since configurations at 
different rows are independent, 



in 



FM\ 



n-1 

n 

1=1 



N 
Mi 



To describe the transitions of this new Markov chain, we need a new 
notion. 

Definition 2.3. Given a configuration C G , the ringing path = 



(Pi , ■ • • , Pn-i) i- s defined as follows: P^li = i and 



>(*) 



P 



(0 



r-1 



P^ + 1, 



i/C 

J 1 

(mod AT) »/ C 



r,P, 



(0 



r,P, 



There is a possible transition between P, 



1 and P^ i+1 ^ /or every row 



r. If we set a = (7_ p (i+i)_ 1 and 6 



C_ p(i+i) , i/ien a 6 —7- 6 a i/ and onh/ i/ 



= o and b = *. This is called a ringing path transition at site i. 

We illustrate this idea with an example in fifrls 2 ^ n Fig ure [2] 

We are now in a position to define the Markov chain on these queues. 
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00*000*0 00*000*0 



00*00 



o * ► 00*0*00* 



*0**G)*** *0**0*** 

Figure 2. The clock rings at row 3 at site 5, causing the 
red ringing path transition P^ = (6,6,5). Only the * at 
row 2 moves. If the clock rang at site 4 instead, then the 
blue ringing path P^ = (5,4,4) would cause no transition. 

Definition 2.4. The Ferrari-Martin multiline queue is the Markov 
chain on SI™ 1 where the dynamics occurs by ringing path transitions that 
take place via exponential clocks of rate 1 at all sites on row n — 1. 

Ferrari and Martin proved a number of remarkable properties of these 
transitions [8]. The first observation is that these transitions are invertible. 
Since they occur with rate 1, this implies that the stationary distribution 
is uniform for all multiline queues. They also showed that each row of the 
multiline queue performs a simple multispecies exclusion process with n = 2, 
where *'s behave like particles and o's behave like vacancies. 

2.3. Bully path projection. There is a formal notion of projection of 
Markov chains called lumpability Lemma 2.5], which we state for 
completeness. Let f2 be the state space of a Markov chain with transition 
matrix M. If f2 can be partitioned into equivalence classes, denoted [a], so 
that M(x, [y]) = M(x', [y]) whenever x' 6 [x], then the resulting process 
on the equivalence classes is also a Markov chain. Then £1 is said to be 
lumpable with respect to the equivalence relation. 

Ferrari and Martin [8] proved that there is a way to lump the Markov 
chain onto Q m in a way so that one obtains the multispecies exclusion 
process. We explain this procedure now. This projections has the property 
that if one looks only at the first j rows of the multiline queue, then one 
obtains the stationary distribution of the multispecies exclusion process with 
m! = (mi, . . . ,rrij,N — Mj ) . 

We now define this projection B : ^l^ M — > fi m . It is defined recursively 
by paths that we call bully paths. A bully path is a path going through 
the multiline queue which always moves rightwards or downwards and which 
contains exactly one * from each row. Moreover, it moves downwards from 
a given row j if and only if it has encountered a * in row j that has not 
already been part of another bully path. 

We start by defining bully paths starting at each of the *'s in the first 
row. The order in which we run these paths does not matter. We mark by a 
1 all the mi columns where the bully paths end at the bottom row. Now, all 
the *'s in the second row which were not part of the bully paths of the first 
class particles are now our second class particles and we run new bully paths 
beginning at each of them. We then mark all the 7712 columns in the bottom 
row that these paths reach by a 2. We continue in this way for all the *'s in 
all the rows. The m,- bully paths starting in row j gives a j on row n. At the 
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end, we have a configuration on row n of entries from {1, . . . , n}. This is the 
required configuration in Q m . We call this projection B. We remark that 
B is well-defined in the sense that the order of the bully paths starting at a 
given row do not matter. For example, both the multiline queues in Table [2] 
project to 34134231. A more complicated example is given in Figure [3j 



2 



Figure 3. An example of a multiline queue with its set of 
bully paths for N = n = 6. 



We concentrate on following conjecture of Lam and Williams [10], where 
N = n. Let wq denotes the reverse permutation n, n — 1, . . . , 2, 1. 

Conjecture 2.5. (Lam and Williams [10] ) With the normalization v(wo) = 

(VL(V) (1) 



X 



. . . X 



n-2 



every stationary weight v(ir) is a polynomial with pos- 
itive integer coefficients. 

Ferrari and Martin [8] have already shown that the stationary distribution 
on multiline queues is uniform when all the Xj's are equal to 1. This means 
that to find the stationary weight of a given' permutation it suffices to 
count the number multiline queues that project to that permutation. Let 
us point out that there is no known easy way of doing this. In particular, 
the conjecture by Lam [S] about the stationary weight of the identity is not 
proved. It would amount to proving that there are Y\S=i (T ) m ultiline 
queues that project to the identity. 

Our aim is to define a Markov chain on the multiline queues such that 
the stationary weight in the general setting is a monomial, and the sum of 
the monomials corresponding to a certain multipermutation under the pro- 
jection B leads to the stationary distribution of the multispecies exclusion 
process with the transition rules given by §2§. This would give a combinato- 
rial proof of Conjecture |2.5[ and moreover a generalization of that conjecture 
to multipermutations. 

We have not settled their conjecture in full generality, but we do have 
partial results which will be presented in the remaining three sections. We 
now present a conjecture describing the stationary weight of each multiline 
queue. Let a be a configuration of length N with n different classes of 
particles. A vacancy o in a is called i-covered if it is traversed by a bully 
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path from an z-class particle (i.e. starting on row i) , but not traversed by 
any path from an z'-class particle, with i 1 < i. For 1 < i < j < n let Zij(a) 
be the number of o's on row j of a that are i-covered. Let Cj := N — Mj be 
the number of vacancies o on row j, and let Cj := Y17=j+l Ci ^ e ^ ne number 
of vacancies below row j. Our main conjecture is as follows. 

Conjecture 2.6. There is a Markov chain on the state space ofQ™ 1 suc h 
that the stationary weight of any configuration a is given by 

l<i<j<n ^ % ' 

Our results in Section [3] and Section [4] support this conjecture. In those 
sections we use the ringing path transitions of Ferrari-Martin. It is however 
possible to show that one cannot in any simple manner give weights to the 
ringing path transitions when n > 4 in general and get a monomial as weight 
for each multiline queue. 

Example Consider the multiline queue a in Figure [3j Since it is a per- 
mutation we have Cj = P 1 ^)) for all j. The only non-zero Zjj's are 



£1,3 = 2, 2^4 = zi,5 = 22,3 = 1- Conjecture 2.6 states that 



/ \ 10 6 3 X3 \ X A \ X 5 \ X 3 \ 6 5 6 2 
v{a) = Xl W 4 [-) {-){-){-)= ^1 W4«B. 

We have checked these conjectures for all possible configurations up to 
size N = 6 on Maple™. We did these by comparing two quantities. First, 
we calculated the stationary weight for each multipermutation by looking at 
the inhomogeneous multispecies exclusion process. Secondly, we calculated 
the weight for each configuration by adding the stationary weight as given 



by Conjecture 2.6 of each multiline queue which projected to it. 



3. Three species Ferrari-Martin model 



Our first result is a proof of the natural generalization of Conjecture 2.5 
for multipermutations with three classes, i.e. n = 3. To state the result, 
we will simplify notation in the following way. Let a be a configuration in 
^mfm 2 m 3 wri ose projection Ba = tt with 7Tj = 3. We say that the 3 is 
covered if the o on row 2 and column i is 1-covered, i.e. if there is a bully 
path passing through that o. If no such bully path exists, we say that the 3 
is non-covered. 

We will use the same transitions as in the Ferrari-Martin model, the 
ringing path transition, see Section [2j which we will assign rates in a special 
way. Since we want our Markov chain to project to the inhomogeneous one 
on f2 m we set the transition rate between two multiline queues projecting 
to different permutations to the same rate as in Q m . This will guarantee 
that the lumpability criterion is satisfied. Thus we define the rate of the 
transition in O,^ 1 from a at column i to take place with rate XB ai when 
Bcti = 1,2. Furthermore, we make the following choice. If Ban = 3, then 
set the rate to be x\ (resp. x^) if the 3 is covered (resp. non-covered). 

The graph of the Markov chain for the multiline queue m = (1, 1, 1) is 
given in Figure |4j 
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Figure 4. All transitions and the corresponding rates for 
the inhomogeneous Ferrari Martin multiline queue with m = 
(1,1,1). The stationary weights are given in parenthesis and 
in red. 

We now collect some observations about these transitions in the following 
lemma. Each of these can be verified without too much difficulty. We will 
write Bon for B(a)i. 

Lemma 3.1. (1) If Ban = 2, the multiline queue at site i is fixed to be 
of the form 

o 
* 

independent of the label at sites before or after i . 
(2) // Baii~± = 1 and Bcti = 2, the Ferrari- Martin transition at site i 
does nothing. 

o/* o 
* * 



1 2 

Independent of whether ax,i—i = ° or *, nothing changes. 

(3) Suppose Bcti, Bcti+i, . . . , Bctj = 3 and Ba%-i, Bctj+i ^ 3. Then 
there is at most one non-covered 3 in this block which can cause a 
nontrivial transition. 

(4) The number of covered 3's can increase (resp. decrease) only with 
transitions starting at sites with non-covered 3's (resp. 1's). 
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Theorem 3.2. Let m = (mi, 1712,1713). With transition rates as above the 
stationary weight of a configuration a £ f^ M is x™ 3_ x\, where k is the 
number of covered 3's in a. This projects to the required stationary distri- 
bution of Ba £ f2 m . 

Proof. We will prove the result by showing that the formula for the sta- 
tionary weight satisfies the master equation at each state. Equivalently, we 
want to show that the total probability leaving a given configuration a is 
equal to that entering it. This is easiest if we represent the corresponding 
multipermutation in block form, 

Ba = i a i2 bl 3 Cl l a2 2 b2 3 C2 • • • r-2 h . 

where we want this representation to be maximal and unique, so that one 
of ai , bi or a can be zero but not two consecutive one of these can be zero. 
Further, suppose that the number of non-covered 3's in a is k, so that the 
stationary weight of a is x™' i ~ k x\. 

Let us first analyze the X2 transitions. These occur if the Ferrari-Martin 
transition starts at a site i such that Bon = 2 or a non-covered 3. By 
Lemma |3.1[ 1), the transition at 2 can only occur if Ba{-\ = 3. Therefore, 
all such transitions occur either inside a block of 3's or at its boundary. 
Let us first look at outgoing transitions from a focussing our attention on a 
block of 3's. If the leftmost 3 in the block is covered, then so are all the other 
3's. This also forces the site to the right of that block to be a 1. Therefore, 
no X2 transition can take place. The only X2 transitions that occur happen 
when the leftmost 3 is not covered. We now have two possibilities. 

(a) There is a transition by a 2, which happens only when all the sites 
in the first row of the block are o's. 

00... |o o| 00... I o o I 

o o ... |o *| y 00... J* o| 

3 3 ... 3 2~ 3 3 ... 2 3~~ 

(b) There is a transition by a 3 within the block. By Lemma |3.1[ 3), 
only one such transition can occur. 

o ... o |o *| ... o o/* o ... o I* o| ... o o/* 

o |o o| o ... o * — y o |o o| o ... o * 

3 ... 3 3 3 ~. 3 T~~ 3 . . . 3 3 3 ~. 3 T 

The important fact is that the (a) and (b) transitions are mutually exclusive 
as well as exhaustive. Therefore, there is one X2 transition for every block 
of 3's, the leftmost of which is not covered and no other X2 transitions. 

Now, let us look at incoming transitions. We want to show that the only 
transitions with an "effective" rate of X2 are exactly one per block of 3's 
with a non-covered leftmost 3. By effective rate X2, we mean contributing 
X2V (a) to the incoming side of the master equation, which happens in tow 
different ways. Either the configuration leading to a has k covered 3's and 
the transition has rate X2 or it has k + 1 covered 3's, i.e. stationary weight 
x ™s-k-i x k+i an( j ^ e transition has rate x\. This time, the possibilities 
depend on the site to the left of the block in a. This is because a nontrivial 
transition at a site occupied by a non-covered 3 comes from a state with 
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fewer covered 3's and a transition at a covered 3 has rate x± and comes 
from a state with the same weight. To differentiate between incoming and 
outgoing transitions, we always place a on the left. 

(A) If the left is occupied by a 2, then the transition happens with rate 
X2 and the number of covered 3's is unchanged. 



2 3 3 ... 3 3 2 3 ... 3 

(B) If the left is occupied by a 1, then the transition happens with rate 
x\ and the number of covered 3's decreases by 1 (since we assumed 
the leftmost 3 of the block of a to be noncovered). 

| ° /* °| • ••• • 1°/* °l ' 

* o| o . . . o i |o *| o . . . o 



1 3 3 ... 3 3 1 3 ... 3 

There is a subtle point here, which we should emphasize. Note that 
the configuration on the right has a o at the first row in the column 
corresponding to 1. If there is a * at that position, that configuration 
may also go to a with rate x\. However, the number of covered 3's 
will not change, and the transition will have an effective rate of x\. 

We have thus shown that there is exactly one incoming transition with 
an effective rate of X2 for every block with a leading non-covered 3. We now 
show that all other incoming transitions have effective rate x%. 

If a non-trivial transition starts at a covered 3 then the effective rate is 
x\ by the discussion before (A) above. If a 1 has jumped left over a 2 or a 
non-covered 3 the effective rate is clearly x\. The last possibility to consider 
is when a 1 has jumped over the 3, and the 3 is covered after the jump, as 
shown below. 

|o/* °/*| ■ ••• ■ 1°/* °/*l ■ 

* o| o . . . o i |o *| o . . . o 



1 3 3 ... 3 3 I 3 ... 3 

This happens with rate x\ and the 3 must have been covered both before 
and after the jump so the effective rate is also x\. Note that the 3 may be 
covered by a * directly above it or somewhere to the left outside this picture. 

To complete the proof, recall that the total number of incoming transitions 
equals that of the outgoing transitions [HE]. □ 



4. Ferrari-Martin model with one first class particle 

Our next result is for arbitrary n with the condition m\ = I. We also 
fix the rates so that we only have one free parameter. Recall that Zij(a) is 
the number of o's on row j of a that are i-covered and that C\ is the total 
number of o's on row 1. Let z\{a) = J2j z i,j( a )- 

Theorem 4.1. In with m\ = 1 and Xi = 1 for 2 < i < n, set the rate 

of the Ferrari- Martin transition from a at column i to take place with rate 
xi when Boti = 1. If Ban = 2, . . . ,n, then set the rate to be 1. Then the 
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stationary distribution of a is v(a) = , where Z m is the normalizing 

•partition function. 

Proof. Let the probability of configuration a be given by v{a). Since it is 
the unique solution of the master equation, it suffices to show that v{a) = 
x l Zl '^ a ) satisfies it. 

Ferrari and Martin [8] showed that there are as many incoming transi- 
tions into any configuration as outgoing ones, say k. Therefore, we need to 
compare the number of incoming transitions with an effective rate of x\ with 
the number of outgoing transitions with rate x\. Since there is only one first 
class particle, the only time we get an outgoing transition with rate x\ is if 
the Ferrari-Martin transition occurs at the unique site i so that Bcti = f. 
There is thus a contribution of [x\ + k — l)v(a) to the outgoing transitions 
from a. 

Now, let us look at incoming transitions to a in the master equation. 
The key observation is illustrated by Figure [5j which is a cartoon for a 
configuration a' along with a generic Ferrari-Martin transition at j to a. 
The only ringing transitions that affect the bully path are at the locations 



j 



Figure 5. The ringing path starting at j intersecting with 
the first class bully path. This transition has rate 1. 



where both paths cross. They cross either in a single vacancy o or in a 
vertical tower of *'s as shown in Figure [5} In the first case the transition 
does not change the bully path. 

When the transition takes place in the latter case, all the *'s at the in- 
tersection will move left, if possible. If any particular * did not move, that 
means there is already a * in the position to its left. Therefore, the vertical 
part of the bully path will shift by one to the left in a. This means that one 
less o is covered at the top of the vertical part and one more o is covered at 
the bottom by the first class bully path. Note that the ringing path could 
intersect the bully path multiple times because these multiline queues live 
on a cylinder. But the argument given above holds true for every vertical 
segment of the intersection of the two separately. Therefore, z\{a') = z\(a), 
and this transition contributes 1 • v{a') = v(a) by assumption. 

We have shown that a transition whose ringing path intersects the first 
class bully path somewhere in the middle does not change Z\ after the tran- 
sition. The only possible way for the transition to affect z\ is if the ringing 
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path intersects the bully path either at the top or the bottom or both. As 
we saw above, it could intersect the bully path at other places, but they 
don't matter. We shall consider these three cases separately. 

(1) Bottom: In this case, the ringing path in a' starts at the site con- 
taining the *, and the transition has rate x\. By the argument 
above, a' has one more covered o than a. Therefore this transition 
contributes x\ ■ v(a') = v(a) and has an effective rate of 1. 

(2) Top: The transition has rate 1 since the ringing path in a' starts at 
a site not containing a 1. But this time, there is one less covered o 
in a' than a and therefore the contribution is 1 • v(a') = Xiv(a). 

(3) Top and Bottom: The transition has rate x\ since the ringing path 
in a 1 starts at the site containing 1. But this time, there is one less 
covered o at the top and one more at the bottom in a' compared to 
q and therefore the contribution is x\ ■ v(a') = x\v (a). 

This shows that the only incoming transitions which have an effective 
rate of x\ are those which intersect the first class bully path at the top. By 
the Ferrari-Martin inverse construction, there is exactly one such transition. 
Thus, there are k — 1 incoming transitions with effective rate 1 and exactly 
one with rate x\. □ 



In this case we may use the theorem to compute the partition function Z m , 
i.e. the divisor for all stationary distributions to get a probability measure 
with the normalization we are using. Recall that we have Mj = Yl)=i m j- 

If we let [k] q := and denote by [k]q® the d'th derivative of [k] q with 

respect to q. Then [k] { q d) = d\ Eto^' OtV" 

Proposition 4.2. In f2„ M with m\ = 1, x± = a and xi = 1 for 2 < i < n, 

the partition function is 

n-l ^(Mj-l) n-l 

Z ™ = N II (M.-Dl =N U hn-jiha^^a), 
j -' { 3 >' i=2 ) 

where is the complete homogeneous symmetric polynomial. 

Proof. The latter equality is an immediate consequence of the definition. We 
will prove the former inequality combinatorially by looking at the multiline 
queues. 

Since we have assumed that mi = 1 and Xi = 1 for i > 2 we need to 
focus on the bully path of the 1. For 2 < j < n — 1, let ej = z±j. For a 
fixed set of < ej < N — Mj for 2 < j < n — 1, such a bully path can be 
completed to a multiline queue by choosing, for each row j, the remaining 
Mj — 1 positions of particles among the N — ej — 1 possible positions. Since 
the 1 has highest priority in every queue this can be done independently for 

every row. Let k := Y17=i e * anc ^ := Yln FM a * ' ■ ^ ne e i' s determine the 
bully path except for the starting position which gives a factor N. For every 
vacancy o that the bully path goes through, the stationary weight gets a 
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factor a 1 by Theorem |4.l[ Together we get 
2®=N-a« Y ( N -^- l \a^...( N - e ^- 1 

0<ej<N-Mj 
k=J2 e j 

To get the partition function, we use Z m = Zm . We get by standard 
manipulations the following formula. The easiest way to see the identity is 
to expand the product below, which gives all the possible terms above. 

*-^n((*: 1 1 ) + (i:V + - + (^:0^) 

Using that C\ = Y27=2 ^ ~ ^i, we can rewrite the formula above as 



a 



j=2 i=0 j=2 i=0 

which can be seen to be the same as the desired formula. □ 
5. A minimal Markov chain for n = 3 



As mentioned after Conjecture 2.6 for re > 4 we do not think that there is 



a simple rule for transition rates in general for the ringing path transitions in 
Q m M that would give each multiline queue a monomial as stationary weight. 
For example, when N = n = 4, assuming all transition rates are of the form 
Xi and all stationary weights are monomials, all master equations cannot be 
satisfied simultaneously. 

We have therefore sought other transitions between multiline queues to 
create other Markov chains for which we might be able to find appropriate 
rates. We will in this section present such a Markov chain Q c for the case 



n = 3. It gives the same stationary weights as Theorem 3.2 and Conjecture 



2.6 and is thus an independent solution of the TASEP for n = 3, parallel 
to the work by Angel pQ. For £l c we believe that it might be possible to 
generalize to higher n, but we have not yet been able to do so in general. 

We call Q c the Multiline coupe process. It will be minimal in the 
sense that there will be no transitions between different configurations cor- 
responding to the same permutation. 

First we will need some definitions. In this section m = (mi, 1112,1713). We 
will divide a word w G £l m into pieces. First we cut w to the right of every 
2 not followed by a 2. Then we will cut each such piece to the right of every 
1 that is not followed by a 1. Every subword thus obtained will start with a 
number (possibly zero) of 3's followed by some consecutive l's or consecutive 
2's. Each such subword will be called a coupe. There are thus four kinds of 
coupes, as can be seen in Figure [6j Each coupe contains exactly one among 
l's or 2's and are naturally called first-class and second-class coupes. We 
call a first-class or second-class coupe with no 3 full. 

In a coupe we will call the leftmost 1 (or 2) the front seat and the 
rightmost letter in each coupe is called the back seat. It is easy to see 
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••• | 3 3 2 | 3 3 1 1 | 3 2 2 2 | 1 | 2 | 3 3 1 1 1 | 2 | ••• 

Figure 6. A decomposition of a circular word into coupes, 
seperated by vertical lines. 

that the only particle in each coupe that can jump is the front seat particle. 
Similarly, the only way to get to a given configuration is if the back seat just 
jumped. The front seat particle can always jump except when the coupe 
consists of only 2's. In that case, the back seat of the coupe to its left is a 
1, and so the 2 cannot jump. Note that a 2 in a front or back seat must 
always be vacant. 

The partition of a word into coupes induces a partition (via the bully 
path projection B) of each configuration corresponding to that word. If in 
a configuration a, a back or front seat has a * in row 1, then it is said to be 
occupied, if it has a o in row 1, then it is said to be vacant. 

Definition 5.1. The transitions in Q c are given by the following rules. 

• Regular jump - If a 1 corresponding to an occupied front seat jumps 
then both the * 's in row 1 and row 2 jumps to the left if possible. 

• Pulling jump - If a 1 or a 2 corresponding to a vacant front seat 
jumps then all * 's in the first row and to the right of the jumping * 
but in the same coupe also moves one step to the left. If the jumping 
1 or 2 corresponding to a vacant front seat is also a back seat then 
all the * 's in the coupe to the right are moved one step to the left. 

Here are two examples of regular jumps: 
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followed by two examples of pulling jumps 
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3. 



Note that a jump might increase the number of coupes by one, decrease 
it by one, or leave it unchanged. Note also that after the pulling jump the 
coupe to the right will have a vacant back seat. 

The rate when a 1 (resp. 2) jumps is x\ (resp. X2) since we want the 
Multiline coupe process to be lumpable with respect to the bully path pro- 
jection. As an example of the multiline coupe process we give 0^ when 
m = (1, 1, 1) in Figure [7J 
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o o * 
o * * 




( Xl ) * * o * * o( X2 ) 

12 3 1 2 3 

Figure 7. All transitions and the corresponding rates for 
the inhomogeneous Multiline coupe process with m = (1, 1, 
1). The stationary weights are given in parenthesis and in 
red. 

Remark 5.2. The jump rules above could be thought of a sequence of 
Ferrari- Martin jumps, but not in a straightforward way. Maybe proofs or 
generalizations will be easier if viewed that way. Note also that row 1 will 
not behave as a TASEP on its own. It would be interesting to understand 
this process. 

The following lemma describes the reverse, i.e. the incoming transitions 
to a given state. The verification of the Lemma is not difficult and we omit 
the details. 

Lemma 5.3. Let a be a configuration. Transitions leading to a in are 
determined by the back seats in the following way. 

(i) An occupied back seat can always have been the last 1-class particle 
to jump by a regular jump. 

(ii) Any back seat where the coupe to the right is not a full second-class 
coupe and has a vacant back seat could have made a transition via a 
pulling jump. 

First we show that we the coupe Markov chain is irreducible, i.e. that we 
can reach any state from any other state. 

Lemma 5.4. For any m = (mi, 7712,711.3) the Multiline coupe Markov chain 
OS! is irreducible. 
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Proof. Let wq be the word starting with all the 3's, then all 2's and then all 
l's. First we prove that we can get from any state a to every cyclic shift of 
wo. Begin by moving all the l's left until they are all directly to the right of 
a 2. Then every 1 will have *'s in both rows above it. Next, jump with the 
2 and then all l's to get to the next 2 and so on until all 2's are together. 
Then again jump with all the 2's and then all l's performing a rotation until 
the desired cyclic shift of wo is obtained. 

Conversely, we will show that we can get from some cyclic shift of wo to 
any state, or equivalently, that we can perform inverse jumps from any state 
a as described by Lemma [573] , and reach some cyclic shift of wq. 

We begin by locating a first-class coupe X to the left of a second-class 
coupe Y. If X has a vacant back seat we do inverse jumps (ii) with the 
back seat in the coupe to the left of X until the back seat of X is occupied. 
Then do regular inverse jumps (i) with the back seat of X to the right of 
all the 2's in Y. Now, if there are a still a first class coupe to the left of Y 
we repeat this procedure, moving a new 1 to the right of Y. If there are a 
second class coupe to the left of Y we merge the two second class coupes 
into one by inverse jumps of type (ii). Continuing this process will give us 
some cyclic shift of wq. □ 

Theorem 5.5. For any m = (mi, m^^mz), the jumping rules defined above 
for ms a Markov chain fi^ that has the stationary distribution of Theo- 
i.e v(a) = x™ z ~ x\, where k is the number of covered 3's. More- 



3.2 



rem 

over, it projects to a solution of tt r 

Proof. Let a be a configuration of with stationary weight v(a) and which 
is divided into c coupes. We have to check the master equation for a, i.e. 
that the inflow is equal to the outflow at a. We will now show that the 
number of incoming transitions equals the number of outgoing ones. For 
this purpose, it will be convenient to distinguish the coupes. Let c\ (resp. 
C2) be the number of first-class (resp. second-class) coupes. Similarly, let f± 
(resp. f'2) be the number of full first-class (resp. full second-class) coupes, 
and ei = c% — fi. Further, let o\ be the number of first-class coupes whose 
back seat is occupied, and v\ (resp. ^2) be the number of first-class (resp. 
second-class) coupes whose back seat is vacant. Note that C2 = t>2- 

In each first-class coupe exactly one particle can jump, namely the front 
seat. In the second-class coupes the front seat can jump if and only if it is 
not full. Thus, the total number of transitions going out of a is equal to 
ci + e2- By Lemma 5.3 (i) there is exactly one transition coming in for each 



occupied back seat by a regular jump. There are o\ of these. By Lemma 5.3 



(ii), there is also one transition coming in to a for each coupe unless the 
coupe to the right has an occupied back seat or is a full second-class coupe. 
This is given by v\ + v 2 — /2 , and so the number of transitions coming in to 
a is the same as the number of transitions leaving. 

We must now also check the weights. For each second-class coupe there 
is one transition leaving with rate X2 if and only if the second-class coupe 
is not full. Each first-class coupe gives an outgoing transition with rate x\. 
So the outgoing side of the master equation is C1X1 • v(a) + e2X2 • v (a). 

Now to the incoming side of the master equation. Any incoming transition 
corresponding to a regular jump has weight a and the number of 1-covered 
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o's on row 2 has not changed so it will come from a configuration with the 
same stationary weight. The total incoming contribution from regular jumps 
is thus o\X\ ■ v(a). For the incoming pulling jumps there are several cases 
to consider: 

I. A back seat in any coupe with a first-class coupe to the right with 
a vacant back seat. By Lemma 5.3 (ii) this corresponds to an incoming 
transition with weight x\. Because of the definition of a pulling jump, the 
number of 1-covered o's in row 2 do not change (either it is covered both 
before and after, or neither before nor after), so the transition comes from 
a configuration with the same stationary distribution and thus contributes 
x\v{oi). The total contribution is v\x\ ■ v(a). 

II. A back seat in a first-class coupe at position i with a second-class 
coupe to the right that is not full. By Lemma 5.3 (ii) this corresponds to an 
incoming edge with weight x\. This time the number of 1-covered o's in row 
2 does change. Namely the o in row 2 at position i + 1 is not 1-covered in a, 
but was 1-covered in the previous configuration, which thus had stationary 
weight v(a) • The contribution for this edge is thus x\ ■ v(a)^r = X2v(a). 

III. A back seat in a second-class coupe that has another (not full) second- 
class coupe to the right. This incoming transition will have weight xi and 
come from a configuration with stationary distribution v (a) since the num- 
ber of 1-covered o's has not changed. Hence contribution X2v(a). 



The total contribution from cases II and III is C2 — /2- By Lemma 5.3 



these are all possibilities, and it is easy to check that the master equation is 
satisfied for a. 

Finally, we must prove that B lumps 0^ to a solution of Q, m - To see this 
it suffices to note that for every a E with B(a) = w where w has a 
transition to some other word w' , there is a unique transition from a to /3 
in il^j such that B((3) = w' and this is clearly true by the definition of the 
jump rules. □ 
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